PHOTON EMISSION FROM BARE QUARK STARS 



B.G. Zakharov 

L.D. Landau Institute for Theoretical Physics, GSP-1, 117940, 
O 

Kosygina Str. 2, 117334 Moscow, Russia 
< 

^. 

We investigate the photon emission from the electrosphere of a quark star. It is shown 
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r^ I that at temperatures T ~ 0.1 ^ 1 MeV the dominating mechanism is the bremsstrahlung 

t:j- ■ due to bending of electron trajectories in the mean Coulomb field of the electrosphere. 

O I The radiated energy for this mechanism is much larger than that for the Bethe-Heitler 
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bremsstrahlung. The energy flux from the mean field bremsstrahlung exceeds the one 
from the tunnel e^e^ pair creation as well. We demonstrate that the LPM suppression 
of the photon emission is negligible. 



1 Introduction 

The hypothesis of quark stars made of a stable strange quark matter (SQM) [1, 2, 3] 
has been attracting much attention for many years. It is possible that quark stars (if 
they exist) may be (at least in the initial hot stage) without a crust of normal matter [4]. 
Contrary to neutron stars the density of SQM for bare quark stars should drop abruptly at 
the scale ~ 1 fm. The SQM in normal phase and in the two-flavor superconducting (2SC) 
phase should also contain electrons (for normal phase the electron chemical potential, /x, is 
about 20 MeV [2, 5]). In contrast to the quark density the electron density drops smoothly 
above the star surface at the scale ~ 10^ fm [2, 5]. For the star surface temperature T '^ fi, 
say T < lO^^K ~ 1 MeV, this "electron atmosphere" (usually called the electrosphere) 
may be viewed as a strongly degenerate relativistic electron gas [2, 5]. 

From the point of view of distinguishing bare quark stars from neutron stars it is 
of great importance to have theoretical predictions for the photon emission from bare 
quark stars. Contrary to neutron stars (or quark stars with a crust of normal matter) 
the photon emission from quark stars made of a stable self-bound SQM may potentially 
exceed the Eddington limit. This fact may be used for detecting a bare quark star. 
However, the SQM itself is a very poor emitter at T ^ cj^ [6, 7] (here u'^ ~ 20 MeV is 
the plasma frequency of the SQM [6]). At such temperatures the photon emission from 
the quark surface is a tunnel process, and the radiation rate turns out to be negligibly 
small as compared to the black body radiation [6]. However, for the electrosphere the 
plasma frequency, cUp, is much smaller than that for the SQM. For this reason the photon 
emission from the electrosphere may potentially dominate the luminosity of a quark star. 



For understanding the prospect of detecting bare quark stars it is highly desirable to have 
quantitative predictions for the photon emission from the electrosphere. This is also of 
interest in the context of the scenario of the gamma-ray repeaters due to reheating of a 
quark star by impact of a massive comet-like object [8], and the dark matter model in 
the form of matter/antimatter SQM nuggets [9]. 

An obvious candidate for the photon emission from the electrosphere is the 
bremsstrahlung from electrons. It may be due to either the electron-electron interac- 
tion (the Bethe-Heitler bremsstrahlung) or interaction of electrons with the mean electric 
field of the electrosphere. One more mechanism is related to the tunnel e^e~ pair creation 
[4, 10]. The point is that the electric field of the electrosphere should be very strong. It 
may be about several tens of the critical field for the tunnel Schwinger pair production 
Ecr = ml/e [11] (we use units c = h = ks = I)- In this scenario the photons appear 
through e+e^ annihilation in the outflowing e^ wind [12]. 

The bremsstrahlung from the electrosphere due to the electron-electron interaction 
has been addressed in [13, 14]. The authors of [13] used the soft photon approximation 
and factorized the e^e^ — ?> e^e~ cross section in the spirit of Low's theorem. In [14] it was 
pointed out that this approximation is inadequate since it neglects the effect of the photon 
energy on the electron Pauli-blocking which should lead to a strong suppression of the 
radiation rate. The authors of [14] have not given a consistent treatment of this problem 
either. To take into account the effect of the minimal photon energy they suggested some 
restrictions on the initial electron momenta imposed by hand. In this way they obtained 
the radiated energy flux from the e~e~ — )■ e~e~7 process which is much smaller than 
that in [13], and than the energy flux from the tunnel e+e^ pair creation [4, 10]. In 



[15] there was an initial attempt to include the effect of the mean Coulomb field of the 
electrosphere on the photon emission. The authors obtained a considerable enhancement 
of the radiation rate. However, similarly to [13] the analysis [15] treats incorrectly the 
Pauli-blocking effect. Note also that in the analyses [14, 15] the photon quasiparticle mass 
was neglected. As we will show this approximation is clearly inadequate since the finite 
photon mass suppresses the radiation rate strongly. 

Thus the theoretical situation with the photon bremsstrahlung from the electrosphere 
is still controversial and uncertain. The main problem here, which was not solved in 
the previous analyses [13, 14, 15], is an accurate accounting for the photon energy in 
the electron Pauli-blocking. In the present paper we address the bremsstrahlung from 
the electrosphere in a way similar to the Arnold-Moore- Yaffe (AMY) [16] approach to 
the collinear photon emission from a hot quark-gluon plasma based on the thermal field 
theory. We use a reformulation of the AMY formalism given in [17]. It is based on the 
light-cone path integral (LCPI) approach [18, 19, 20] (for reviews, see [21, 22]) to the 
in-medium radiation processes. For an infinite homogeneous plasma (with zero mean 
field) the formalism [17] reproduces the AMY results [16]. The LCPI formulation [17] 
has the advantage that it also works for plasmas with nonzero mean field. It allows to 
evaluate the photon emission accounting for bending of the electron trajectories in the 
mean Coulomb potential of the electrosphere. Contrary to very crude and qualitative 
methods of [13, 14, 15] the treatment of the Pauli-blocking effects in [16, 17] has robust 
quantum field theoretical grounds. Of course, our approach is only valid in the regime of 
collinear photon emission when the dominating photon energies exceed several units of 
the photon quasiparticle mass. Numerical calculations show that even at T ~ 0.1 MeV 



the effect of the noncolhnear configurations is relatively small. 

We demonstrate that for the temperatures T ~ 0.1 ^ 1 MeV the radiated energy flux 
from the e~ — )■ e~7 transition in the mean electric field turns out to be much larger 
than that from the Bethe-Heitler bremsstrahlung. It exceeds the energy fiux from the 
tunnel e^e~ pairs as well. Also, we demonstrate that, contrary to conclusion of [13], the 
Landau- Pomeranchuk-Migdal (LPM) suppression [23, 24] of the photon bremsstrahlung is 
negligible. Our results show that the photon emission from the electrosphere may be of the 
same order as the black body radiation. For this reason the situation with distinguishing 
a bare quark star made of the SQM in normal (or 2SC) phase from a neutron star using 
the luminosity [4, 25] may be more optimistic than in the scenario with the tunnel e^e~ 
pair creation [4]. 

In a short version the results of this work were presented in [26]. In this paper we 
present our results in a more detailed form. The plan of the paper is as follows. In Sec. 2 
we review the basic formulas and approximations. In Sec. 3 we discuss the evaluation of 
photon emission from a given electron in the electromagnetic field of the electrosphere 
which includes both the mean Coulomb field and the ordinary fluctuation field generated 
by neighboring electrons. In Sec. 4 we present numerical results for the radiated energy 
fiux. Sec. 5 is devoted to conclusions. 

2 Basic formulas and approximations 

As in Refs. [4, 13, 14] we use for the electrosphere the model of a relativistic strongly 
degenerate electron gas in the Thomas-Fermi approximation. In this approximation the 



local electron number density reads ne{h) = fi^ (h) / dn'^ , where h is the distance from the 
quark surface. The h dependence of the chemical potential is governed by the Poisson 
equation for the electrostatic potential V = fi/e. For h > this gives [2, 5] 
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where H = ^37r/2a/fi{0), a = c^/Att. 

We assume that the electrosphere is optically thin. This means that the photon 
absorption and stimulated emission can be neglected. In this regime the luminosity may 
be expressed in terms of the energy radiated spontaneously per unit time and volume, 
Q, usually called the emissitivity. In the formalism [17] the emissitivity per unit photon 
energy w at a given h can be written as 

dQ{h,u)_u{k)dk /^P^^(^)[i_,^(^.)]g(^_^)^^(P>^)^ (2) 



dcu 4:7[^ du J p dxdL 

where k denotes the photon momentum, E and E' are the electron energies before and 
after photon emission, np{E) = [exp{{E — fJi)/T) + 1]^^ is the local electron Fermi dis- 
tribution (we omit the argument h in the functions on the right-hand side of (2)), and 
X = k/p is the photon longitudinal (along the initial electron momentum p) fractional 
momentum. The function dP/dxdL in (2) is the probability of the photon emission per 
unit X and length from an electron in the potential generated by other electrons which 
includes both the smooth collective Coulomb field and the usual fluctuating plasma part 
related to the field generated by the neighboring electrons. Note that the formula (2) 
accounts for photons emitted to all directions, because in the optically thin electrosphere 
practically all the photons radiated to the hemisphere directed to the quark surface will 



be reflected either in the electrosphere (at the level with oo^ = u) or from the quark sur- 
face. Only the photons with cj > w^ ~ 20 MeV may be absorbed in the quark matter. 
However, such photons are not important at T < 1 MeV considered in the present paper. 
For the above reasons, it would be incorrect to exclude the photons emitted towards the 
star surface, as was done in [14]. 

Our basic formula (2) assumes that the photon emission is a local process, i.e. the 
photon formation length ^ (we denote it //) is small compared to the thickness of the 
electrosphere. Evidently, only in this case one can define a local emissitivity. Note that 
Eq. (2) defines the rate of photon production at a given photon energy, which remains 
constant during the photon propagation in the electrosphere. The photon momentum in 
this process changes adiabatically according to the photon quasiparticle dispersion relation 
in the electron plasma. Also, the formula (2) assumes that on the scale ~ If the electron 
trajectories are smooth. It means that besides the evident condition If <^ Rm {Rm is 
the radius of curvature of the electron trajectory in the mean field) the typical scattering 
angle related to the random walk of an electron due to electron-electron interaction should 
be small as well. One can show that these conditions are satisfied for the electrosphere. 
An important consequence of the smoothness of electron trajectories at the scale ~ If is 
the longitudinal factorization of the Pauli-blocking factor [1 — nplE')] for the final state 
of the radiating electron in (2). Namely the fact that the trajectories are smooth in the 



^Physically, the photon formation length (sometimes called the coherence length) is a longitudinal 
scale at which the photon and electron wave packets become separated. It appears naturally in the 
LCPI approach [18, 21] formulated in the coordinate space as a dominating scale of the integrals in the 
longitudinal coordinate. 



process of the photon emission allows one to neglect the statistics effects in treating the 
small angle scattering. Indeed, the typical space scale for the soft fluctuating modes of 
the electromagnetic fleld is about the inverse Debye mass I/itid ~ l/e/x. This scale is 
much larger than the typical separation between electrons ~ l//i. For this reason from 
the point of view of the electrons with energy ~ /i the soft electromagnetic fleld at the 
space scale ~ I/ttld ^ 1/yU can be viewed as a uniform fleld at the scale ~ l//i. In a 
uniform fleld all electrons in the same spin state scatter in the same way, and the small 
angle scattering leads simply to some shift of the distribution function in the momentum 
space. Any statistics effects will be suppressed by some power of electron charge e. 
Calculations within the real time thermal fleld theory performed in [16] corroborate this 
physical picture of the coUinear photon emission. 

In our approximation of optically thin medium the differential radiated energy flux 
from the electrosphere, dF/du, in terms of the emissitivity reads 

dF^f'-^j^dQih,u)_ (3) 

doo Jo doo 

For the chemical potential (1) the /i-integration in (3) can be approximated by the inte- 
gration over fi as 

dF _ ^ /-^(o) dfidQ{h{fi),u) 



doo V 2a Jfimin /U^ doo 

with fimin = jJ'ihmax)- In numerical calculations we take Umin = ^rrie. Of course, the 
relativistic approximation we made is not good at // ~ rrie. However, the contribution of 
this region is small, and the corresponding errors are not big. 



3 Calculation of dP/dxdL 

The essential ingredient of Eq. (2) is the probabihty distribution dP/dxdL for the photon 
emission in the electromagnetic field of the electrosphere. Due to presence of the product 
np{E)[l — np{E')\ in (2) the emissitivity is dominated by the photon emission from the 
electrons near the Fermi surface with p ~ /i ^ nif.. This allows one to use for the photon 
spectrum dP/dxdL the quasiclassical relativistic formulas. In this work we evaluate this 
spectrum within the LCPI formalism [18, 21]. In this approach it can be written as 

dP 7 

^-^ = 2ReJ dCgix) [/C(p2,e|Pi,0) - /C,(p2,e|Pi,0)] \p^^p^^,- (5) 



Here g is the spin vertex operator given by 

with gi{x) = a(l — x + x^/2)/x and g2{x) = amlx^/2M^{x), M{x) = px{l — x), /C is the 
Green's function for a two-dimensional Schrodinger equation with the Hamiltonian 

where Lq = 2M{x)/e'^, e^ = m^x"^ + (1 — x)'m'i {m^ is the photon quasiparticle mass), 
the form of the potential v will be given below. In (5)- (7) p is the coordinate transverse 
to the electron momentum p, the longitudinal (along p) coordinate ^ plays the role of 
time. The K^ in (5) is the free Green's function for f = 0. Note that at low density and 
vanishing mean field the quantity Lq coincides with the real photon formation length Ij 
[18] which characterizes the dominating scale in the .^-integration on the right-hand side 
of (5). 
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The potential in the Hamiltonian (7) can be written a,s v = Vm + v f . The terms Vm 
and Vf correspond to the mean and fluctuating components of the vector potential of the 
electron gas. Note that when If is small compared to the scale of variation of /x (along the 
electron momentum) one can neglect the ^-dependence of the potential v in evaluating 
dP/dxdL. The mean field component is purely real Vm = ~xfp with f = edV/dp [21, 27]. 
It is related to the transverse force from the mean field. Note that, similarly to the 
classical radiation [28] , the effect of the longitudinal force along the electron momentum 
p is suppressed by a factor ~ (rrie/E)^, and can be safely neglected. The term Vf can be 
evaluated similarly to the case of the quark-gluon plasma discussed in [17]. This part is 
purely imaginary f/(p) = —iP{xp), where 



P{p) = e' J d^[G{C,0^,O-G{^,p,O], 



G{x — y) = u^UiyD^'^, D^^" = {A'^ (x) A"^ (y)) is the correlation function of the electro- 
magnetic potential (the mean field is assumed to be subtracted) in the electron plasma, 
u^ = (1, 0, 0, —1) is the light-cone 4-vector along the electron momentum. Note that the 
function P{p) is gauge invariant by construction, and one can use D^'^ in any gauge. The 
formula (8) may be rewritten as (below we replace the argument of P{p) by p = \p\ since 
P{p) does not depend on the direction of the vector p) 

2 

Pip) = J^ I rfq±[l - exp(^q^p)]D(q^) , (9) 

where the function D in terms of the correlator G in momentum representation reads 

oo 

D{q±) = 7^ dqodqJ{qo - qz)G{qo, q^, q^) . (10) 

ZTT J 
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The function -D(q^) may be expressed in terms of the longitudinal and transverse pho- 
ton self energies IIl,t- We use for them the formulas of the hard dense loop approximation 
(HDL) [30, 31]. The details of the calculations are given in the Appendix A. 

The function P{p) has been introduced, for the first time, in the context of the prob- 
lem of propagation of relativistic positroniums through amorphous media [29], where the 
atomic size plays the role of the inverse Debye mass. In our approach the function P{p) 
contains all information about the electron-electron interaction which is necessary for de- 
scription of multiple scattering of a given electron in the fluctuating electromagnetic field 
generated by other electrons. In particular, all the Pauli-blocking effects in the process 
of electron multiple scattering are automatically accumulated in P{p). It is worth noting 
that in the approximation of static Debye screened scattering centers the function P{p) 
reduces to na{p)/2 [17], where n is the number density of the medium, and 

a{p) = 8a' / dq \ ' ^'^ = ^^^^ [^ " P^nK^ipmn)] (11) 

J \fl -r 1 1 I'D) ''''D 

is the well known dipole cross section for scattering of an e^e~ pair of size p on the Debye 
screened scattering center (in (11) Ki is the Bessel function). In the static approximation 
at p ^ l/rrij:, one can obtain from (11) P{p) ~ nCp''/2, where C ~ AT[a''\n{2/pm£,) is a 
smooth function of p. In the limit p ^ I/ttlj:, the function P{p) in the HDL approximation 
also becomes almost quadratic. 

The quadratic approximation P{p) oc p^ in the LCPI approach is equivalent to the 
Fokker-Planck approximation in Migdal's approach [21]. It is not very accurate but rea- 
sonable for bremsstrahlung in ordinary materials. In this case the dominating p-scale is 
~ l/rrieX, and the spectrum is controlled by behavior of P{p) at the scale ~ l/rrie which 
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is much smaller than the screening radius ~ Xjavfif^Z^I^ (here Z is the atomic number). 
In the case of the relativistic electron gas the situation is quite different. Now, in the 
dominating p- region, the argument of Pi^p) is ~ (0.1 ^ 2)/mD- In this region P{p) is 
essentially non-quadratic. It is well seen from Fig. la, in which we plot the results of 
numerical calculations of P{p) for several values of the ratio T/mo. The results are pre- 
sented in a dimensionless form. For comparison in Fig. la we also show the predictions of 



the static approximation at T = (when ttij? = pJAa/ir) obtained with the dipole cross 
section (11). One can see that at p ~ (0.1 -^ 2)/m£, the function P{p) is almost linear in 

P- 

To demonstrate the relative effect of the longitudinal and transverse modes in Figs. lb,c 
we show separately the contributions related to 11^ and 11^. One sees that at p < l/mo 
the L and T contributions are close to each other. However, at p > 2/mu the longitudinal 
part flattens, while the transverse magnetic one continues to grow (for T/m£, not very 
close to zero). This growth of the transverse part is a consequence of the well known 
absence of the static magnetic screening in the electron plasma. Note, however, that from 
the point of view of the photon emission the growth of the magnetic contribution with p 
is not important since the photon spectrum is dominated by p < 1/e ~ l/m^,. 

The growth of P{p) with temperature is due to the presence of the Bose-Einstein 
factor in the function D (A.l). From Fig. la one can see that the prediction of the HDL 
approximation at T ^ Tno-, similarly to the static model, flatten at p > Ijmj). However, 
the static model exceeds the HDL prediction by a factor ~ 2.5. The fact that at T = the 
static approximation overestimates Pi^p^ is quite natural, since the Pauli-blocking effects 
reduce the effective number of the scatterers. Note, however, that it would be incorrect 
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to interpret the growth of P{p) with temperature as an artifact associated only with 
the decrease of the Pauh-blocking at high temperatures. The function P{p) in the HDL 
approximation accumulate all the collective effects in the soft modes of electromagnetic 
field in the electron plasma at the momentum scale ~ m£, <^ fi. In particular, it accounts 
for the temperature dependence of the density of the plasmon excitations. Note that, 
physically, the appearance of P{p) is due to Landau damping of the longitudinal and 
transverse modes. 

It is worth noting that the collective effects cannot be accounted for consistently in 
the naive modification of the photon propagator in the amplitude of elastic e~e~ — ;> e~e~ 
scattering as was assumed in [13]. One of the consequence of the inadequacy of this pre- 
scription is a strong overestimate of the magnetic contribution in [13]. It is connected 
with the 1/6^ {6 is the scattering angle) behavior of the magnetic contribution to elastic 
e^e^ — )■ e^e^ cross section. To perform the ^-integration the authors of [13] introduced 
some minimal momentum transfer. In contrast to [13] the magnetic contribution to the 
function -D(qi^) is oc 1/q^ at q_|_ — )■ 0^, and the q_|_-integration in the formula for P{p) 
(9) converges at small q^. This change in the small angle behavior of the magnetic 
contribution in the prescription of [13] and in our approach is connected with the dynam- 
ical magnetic screening which was not consistently accounted for in [13]. In principle, 
physically, it is evident that the concept of the elastic e^e^ — )> e^e^ amplitude itself is 
ill-defined for the momentum transfer <'mD, where the collective effects come into play. 



^Thc same occurs in the hard thermal loop approximation for a hot relativistic plasma with zero 
chemical potential [32]. Note, however, that a very elegant formula for the analog of our function £'(q_L) 
obtained in [32] is not valid for a strongly degenerate electron plasma. 
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Note that in terms of P{p) the transverse momentum broadening distribution of an 
electron propagating a distance L through the electron gas can be written as 



^(q±) = 7^ J dpexp [iq^p - LP{p)] 



(12) 



This formula looks like the prediction of the eikonal approximation which neglects the 
variation of the electron tranverse coordinate. However, the path integral calculations 



show that it is valid beyond the eikonal approximation as well. 
Let us turn to calculation of the spectrum with the help of (5). Treating Vf as a 
perturbation one can write 

^(6, p2lCl> Pi) = ^m(6, P2I6, Pl)~^ J d^dp]CmiC2, p2\^, p)Vf{p)ICm{C, p\Cl, Pi) + • • • , 



(13) 



where /C^ is the Green's function for Vf = 0. Then (5) can be written as 



dP dPrr, dP, 



f 



dxdL dxdL dxdL 



(14) 



Here the first term on the right-hand side comes from the /Cm — /C^, in (5) after representing 
/C in the form (13). It corresponds to the photon emission in a smooth mean field. The 
second term comes from the series in Vf in (13). This term can be viewed as the radiation 
rate due to electron multiple scattering in the fluctuating field in the presence of a smooth 
external field. 

The analytical expression for the Green's function for the Hamiltonian with a constant 
force is known (see, for example [33]). In our case the /C^ can be written as 



^m(6,p2l6,Pi; 



M 

2^ 



exp < I 



M{p,-p,f xef-(P2 + Pi) xH^e i 



2^ 



24M L 



(15) 
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with ^ = C,2 — ^i- With this expression from (5) after simple calculations one can obtain a 
spectrum similar to the well known quasiclassical synchrotron spectrum [34] which can be 
written in terms of the Airy function Ai{z) = ^J^Ki/-i{2z^^'^ /3) (here i^i/3 is the Bessel 
function). In the case of interest, for a nonzero photon quasiparticle mass it reads [27] 

U/lrn, 



-Ai {n) + b dyAi{y), (16) 

K Jk 



dxdL 

where a = —2e^gi/M, b = Mg2 — e'^gi/M, k = e^ /{M'^x^fy/^. Inspecting the longitudinal 
integrals for the photon radiation in an external field one can find that the effective 
photon formation length for the mean field mechanism is given by L^ ~ iiiin(Lo, Lm), 
where L^ = (24M/a;^f Y^^ [27]. A similar estimate can be obtained from the criterion of 
separation of the photon and electron wave packets. Note that the analytical expression 
for the Green's function for the oscillator with a constant force is known as well (see [33]). 
Making use of this Green's function one can obtain for P{p) oc p^ the radiation rate in 
the form given in [35], where Migdal's approach within the Fokker-Planck approximation 
was generalized to the case with an external field. The formulas of [35] were used in [15]. 
However, as already noted the approximation P{p) oc p^ is clearly not adequate for the 
electrosphere. 

Let us discuss now the fiuctuation component dPf/dxdL. We represent it in the form 

dPf dPf" dPf'''^ ^ ^ 
^ - — ^— ' — ^ (17) 



dxdL dx dx 

where the first term on the right-hand side corresponds to the leading order in expansion 
in Vf in (13), and the second one to the sum of the higher order terms. The dPf^ /dxdL 
is the analog of the Bethe-Heitler spectrum in ordinary materials, while the dPf^^ /dxdL 
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describes the LPM correction. For the Bethe-Heitler term one can obtain from (5), (13) 

-^ = 2 J dpWix,p,i)Pipx) , (18) 



W{x, p, f) = -Re g{x)(^{x, p, p^, f)(^{x, p, p^, f) ^ _ „ _^ , (19) 

$(a;,p,p',f)= r d^lCm{p,0\p\O- (20) 

Note that for a nonzero f the function W cannot be viewed as a probabihty density for 
the I'je) Fock component of the physical photon (it is even not positively defined). This 
is connected with the fact that in an external field the \'~fe) Fock component is not stable, 
and decays through the tunnel transition into free photon and electron. The analog of 
the representation for the LPM correction derived in [19] for nonzero mean field reads 



dPf 



dx ='^^'^9ix) J d^ J dp^{x,p,p2,i)Pipx)^{x,p,p^,i,^)\^ ^^^^, (21) 



where the function $(a;, p, p^, f, ^) is the solution of the two-dimensional Schrodinger 
equation with the Hamiltonian (7). The boundary condition for $(a;, p, p^, f, .^) is 
^{x, p, pi, f, 0) = $(x, p, pi, i)P{px) . 

In the case of zero f the function W may be written as a density for the \'~fe) Fock 
state 

W{x,p) = WMx,pA>^^})\\ (22) 

where \E'(x, p, {Aj}) is the light-cone wave function for the e — )■ •ye' transition, {Aj} = 
(Ae,Ae',A^) a set of helicities. Note that contrary to the case f 7^ now, due to the 
azimuthal symmetry of the Hamiltonian, the light-cone wave functions have definite az- 
imuthal quantum numbers. The LPM correction in this case also can be written in terms 
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of the light-cone wave functions. The results is similar to that for ordinary materials 
[19, 21] 



dPf 



LPM 



ReY^jdijdp ^* (a:, p, {A,})P(px)<l(x, p, {A,}, • (23) 



The boundary condition for <l>(a;, p, {Aj}, ^) is now $(a;, p, {Aj}, 0) = ^{x, p, {Ai})P(pa;) . 
The light-cone wave functions appear in the formulas (22), (23) from the ^-integrals in (5) 
and (13) of the Green's function /C„j and action of the vertex operator written in terms 
of the helicity projectors as was done in [17]. 

The formulas for the light-cone wave functions are given in the Appendix B. Making 
use of the formulas given there one can obtain for the probability distribution W for the 
e — )■ 7e' transition at f = 



a {[1 + {I- X] 
2^ I Ic 



W{x, p) = -^ '±^^^-^e'Kl{pe) + x'mlKl{pe) , (24) 



where K^i are the Bessel functions. Due to exponential decrease of K^i in (24) the 
dominating p scale in the formula (18) for the fluctuation term is ~ 1/e. 

For nonzero f the azimuthal symmetry is absent. This makes the problem considerably 
more complicated. In the present work we have first calculated the spectrum dPf/dxdL 
for f = 0. We observed that the LPM correction in (17) is negligible as compared to 
the Bethe-Heitler term. Also, the Bethe-Heitler term itself turns out to be much smaller 
than the mean field term dPm/dxdL. It is clear that a nonzero f will make dPf/dxdL 
even smaller. For this reason an accurate calculation of the fiuctuation term for nonzero 
f does not make much sense. We have taken into account the effect of the transverse 
force using qualitative arguments based on the estimates of the coherence lengths with 
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and without tranverse force. The mean field should suppress the coherence length. The 
suppression of the radiation rate should be approximately the same [36]. Thus, the mean 
field suppression factor can be written as the ratio of the formation lengths with and 
without the mean field. The coherence length in the presence of the mean field is ~ L^. 
Without the mean field in the regime of weak LPM suppression the coherence length is 
given by Lq. So one has the mean field suppression factor Sm ~ -^m/-^o- Note that due 
to reduction of the effective formation length the LPM effect should become even smaller 
for a nonzero mean field. 

To illustrate the relative contributions of the mean field and fiuctuation mechanisms to 
dP/dxdL we plot them in Fig. 2 for /i = 10 MeV and T = 0.2 and T = 1 MeV. The mean 
field part shown in Fig. 2 corresponds to the spectrum averaged over all directions of the 
electron momentum. The fiuctuation contribution has been calculated without the mean 
field suppression factor. Note that we perform calculations with the /c-dependent photon 



quasiparticle mass extracted from the relation mi = YiT{,\ k"^ + ^, k) . This gives m. 



rising from m£)/-\/3 at A; ^ mn to m/j/v^ at A; ^ m£, with the Debye mass 



m 



D 



^(//^ + 7r^T^/3). From Fig. 2 one sees that the fiuctuation contribution is suppressed 
by a factor ~ 10"^. To illustrate the role of the finite photon quasiparticle mass we 
presented in Fig. 2 the results for zero m^ as well (thin curves) . It is seen that the photon 
mass suppression (called usually the Ter-Mikaelian effect) is very strong at small x. The 
effect is especially dramatic for the fiuctuation part where the well known 1/x form of the 
spectrum is changed into oc x. This effect was ignored in the analyses [14, 15] where the 



■^We ignore the influence of the medium effects on me [37] since the photon bremsstrahlung in the 
region a; <C 1, which dominates the emissitivity, is not very sensitive to the electron quasiparticle mass. 
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massless formulas have been used. The results shown in Fig. 2 indicate clearly that the 
massless approximation is inadequate. 

As previously mentioned, our calculations show that for the fluctuation mechanism 
the LPM suppression is negligible. This is in a contradiction with the analysis [13] where 
the authors find a very strong LPM suppression (about ~ 1/300 at the photon momentum 
k = 0.5 MeV for electron energy 10 MeV). For calculation of the LPM suppression 
the authors of [13] have used Migdal's formulas with zero photon mass putting there 
Z = 1 . However, one can easily show that Migdal's formulas become inadequate for the 
electrosphere. We explain this in the language of the LCPI approach. Migdal's approach 
[24] corresponds in the LCPI formalism to quadratic parametrization P{p) ~ nCp'^/2. As 
described above, this approximation is not accurate for electrosphere, but nevertheless it is 
suitable for our qualitative analysis. In the quadratic approximation the Hamiltonian (7) 



takes the oscillator form with fi = J —inC x^ / M (x) . The LPM suppression factor, Slpm, 
can be written in terms of the dimensionless parameter rj = \Q\Lo [18, 21]. The LPM 
suppression becomes strong at i] ^ I. In this limit Slpm ~ 3/r7v^ [18]. The LPM effect 
is negligible for rj '^ 1 when Slpm{i]) ~ 1 — lQi]'^/21 [18], Note that even at ?] ~ 1 the LPM 
suppression is relatively small since Slpm{^) ~ 0.86. A very strong suppression obtained 
in [13] is mostly due to the neglect of the photon mass. The finite photon mass reduces 
strongly the Lq and correspondingly the parameter t] (about a factor ~ 400 for k = 0.5 
and p ~ 10 MeV). Also, for the electrosphere there is no the well known large Coulomb 
logarithm ln(l/Q;) ~ 5 (which comes from the logarithm in the dipole cross section [20]) 
in the \Q\, which is present in Migdal's formulas derived for ordinary materials. Both 
these effects reduce drastically the value of rj for the electrosphere as compared to that 

19 



in Migdal's approach. As a result, the LPM suppression in the electrosphere turns out to 
be neghgible. 

4 Numerical results and discussion 

In this section we present the numerical results for the emissitivity and radiated energy 
flux. The results were obtained with some modification of the spectrum dP/dxdL in the 
noncollinear region. As we mentioned earlier, the collinear approximation we use becomes 
invalid for very soft photons with k < m^. In this region the formalisms [16, 17, 18] do 
not apply. In particular, the LCPI approach [18], which assumes that the transverse 
momentum integration comes up to infinity, should overestimate the photon spectrum at 
k < m.y. To take into account this effect (at least, qualitatively) in calculating the radiated 
energy fiux we multiplied dP/dxdL by the kinematical suppression factor Skin{k) = 1 — 
exp(— /c^/m^). This factor does not give a big effect. It suppresses the radiated energy 
by ~ 10 - 15% at T ~ 0.1 ^ 0.2 MeV and ~ 1 - 2% at T ~ 1 MeV. This says that the 
errors from the noncollinear configurations are small. 

In Fig. 3 we show the emissitivity for /i = 5 and /i = 10 MeV evaluated for T = 0.2 
and T = 1 MeV as a function of u. One sees that the contribution of the mean field 
emission (thick sohd line) exceeds the fiuctuation emission without mean field suppression 
(dashes) by a factor ~ 10^. The mean field suppression gives additional reduction of the 
fiuctuation contribution (thin solid line) by a factor ~ 3—4. Note that in our quasiclassical 
approximation at a given fi there is no photon emission at u < Up. For this reason the 
differential emissitivity shown in Fig. 3 vanishes abruptly at u = u^ = m^{k = 0). From 
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Fig. 3 one can see that, despite the Pauh-blocking suppression, even at T = 0.2 MeV 
the contribution of energetic photons with energy about several units of Up is important. 



This demonstrates that the restriction on the photon energy u < J(^p^ + Tnl imposed by 
the authors of [13] is clearly inadequate. 

In Fig. 4 we plot the differential radiated energy flux dF/du for /i(0) = 10 and /i(0) = 
20 MeV obtained for T = 0.2 and T = 1 MeV. For the fluctuation contribution we 
show the results with and without the mean field suppression factor Sm- For comparison 
the black body spectrum is also shown. The mean Coulomb field of the electrosphere 
reduces the fluctuation term by a factor ~ 3 — 4. From Fig. 3, 4 one can see that 
the relative contribution of the fluctuation mechanism is very small compared to the 
mean field emission. Thus, in some sense we have a situation similar to that for photon 
radiation from an atom with large Z. Note that the form of the spectrum for the mean 
field mechanism is qualitatively similar to that for the black body radiation. 

In Fig. 5 we show the total energy flux F = J^ dudF/du scaled to the black body 
radiation as a function of temperature. For comparison we also plot the predictions for 
bremsstrahlung obtained in [13, 14, 15]. We also show there the energy flux from the 
e~^e~ pair production [4, 10]. We define it as 

F± = / dhQ^{h) ^ J^ / 4Q±(/^(/i)) . (25) 

Here Q± is the energy flux from e^e~ pairs per unit time and volume. We write it as in 
[4, 10] Q± = Ee+e-dNe+e- /dtdV , where E^+e- ~ 2(me + T) is the typical energy of e"'"e~ 
pairs, and dNe+e-/dtdV the rate of e^e~ pair production per unit time and volume given 
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by 



dtdV 27r3 V TT ^ V T 

with ^ = ^J^, and the function J is defined as in [10] 



,, , x3ln(l + 2/a;) vrV 



-exp -— ^ J(0 (26) 



3(1 + 0.074a;)3 6(13.9 + a;)^ ' 

From Fig. 5 one sees that in the region T ~ 0.1 ^ 1 MeV the mean field photon emission 
exceeds considerably both the fiuctuation bremsstrahlung and the energy fiux from e~^e~ 
pair production. 

Figs. 4, 5 demonstrate that the energy fiux from the mean field photon emission 
may be of the same order of magnitude as the black body radiation. It says that the 
approximation of optically thin electrosphere is not very good, and the photon absorption 
and stimulated emission may be important. However, since the radiation rate we obtained 
does not exceed the black body limit, they can not modify strongly our results. Note that 
the authors of [15] obtained for T < 1 MeV the energy fiux considerably exceeding the 
black body limit. This can be seen from Fig. 5, where the results of [15] at fi{0) = 20 MeV 
are shown. The authors of [15] claim that the electrosphere may radiate stronger than a 
black body. This statement is obviously incorrect. The violation of the black body limit 
in [15] is just a signal that the thin medium approximation becomes inadequate at high 
emissitivity. As far as a very large emissitivity obtained in [15] is concerned, as we already 
mentioned, it may be due to incorrect description of the Pauli-blocking and neglect of the 
photon mass. 

As we mentioned earlier, our assumption that the photon emission is a local pro- 
cess is valid if // ~ Lm ^ L^i, where Lei is the typical scale of variation of the po- 
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tential Vm along the electron trajectory. For the chemical potential (1) it can be de- 
fined as Lei ~ H fi{0) / fi{h) cos 6 , where 6 is the angle between the electron momen- 
tum and the star surface normal. Evidently, the contribution of the configurations with 
Lm ^ Lei into the photon spectrum will be suppressed by the finite-size suppression factor 
Sfs ~ ram{Lei, Lm) / Lm- We have checked numerically that this suppression factor gives 
a negligible effect. This justifies the local approximation. 

According to the simulation of the thermal evolution of young quark stars performed 
in [25] the temperature at the star's surface becomes ~ 0.2 MeV at t ~ 1 s. However, 
in the analysis [25] the mean field bremsstrahlung was not taken into account. In the 
light of our results one can expect that the cooling of the bare quark star's surface should 
go somewhat faster than predicted in [25]. It is worth noting that in the initial stage of 
the quark star evolution the mean field photon emission can only modify the temperature 
near the star surface. The evolution of the star core temperature is driven by the neutrino 
emission [25] since for an extended period of time the neutrino luminosity is much larger 
than the photon (and e^e~) luminosity [25]. Higher luminosity due to the mean field 
bremsstrahlung increases the possibility for detecting bare quark stars. From the point 
of view of the light curves at t > 1 s it would be interesting to investigate the mean field 
bremsstrahlung for T < 0.1 MeV as well. However, at such temperatures the photon 
emission from the nonrelativistic region of the electrosphere may be important, where our 
formulas become inapplicable. As far as the contribution of the relativistic region fi ^ rrie 
is concerned. Extrapolation of the curves shown in Fig. 5 to T < 0.1 MeV allows one to 
expect that the mean field emission will dominate the energy fiux at lower temperatures 
as well. 
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A remark is in order here on the photon distribution seen by a distant observer. For 
obtained values of the energy flux the radiation cannot stream outward freely. The point 
is that near the star surface the thermalization time in the comoving frame for the e^e~7 
wind is negligibly small as compared to the star radius. This follows from estimates of 
the mean free path, A, related to 7 + e^ — )■ 7 + e^ and 7 + 7 ^ e+ + e^ processes. 
The qualitative calculations give A ~ lO^^cm at T ~ 0.1 MeV and A ~ lO^^cm at 
T ~ 1 MeV. For this reason the e+e^7 wind can be described as a hydrodynamical flow. 
The hydrodynamical description is valid up the the freezeout surface, beyond which the 
radiation streams outward almost freely. For an observer at large distance from the star 
the photon spectrum is close to the black body one with a temperature Text = TfrTjr, 
where Tfr is the wind temperature and Tjr the bulk Lorentz factor of the wind at the 
freezeout level [38, 39]. One can show that for a relativistic wind TfrTfr ~ TjFj [38, 39], 
where Tj is the wind temperature after its thermalization and Fj the bulk Lorentz factor 
of the wind near the star surface. For T ~ 0.1 MeV the electron fraction in the e"'"e~7 
wind after thermalization is small. Simple qualitative calculations give in this case TjFj ^ 
T(3fi;r^/16)^/*^, where k = {F + F±)/Fi,i,. As a plausible estimates one can take Ff ~ 3 
and K ~ 1. Then one obtains T^xt ~ 0.85T. For T ~ 1 MeV the electron fraction in 
the wind after thermalization becomes close to that for relativistic plasma. In this case 
TiTi ^ T{3KTf/Uy/^. Taking k ~ 0.4 one obtains T^xt ~ 0.5T. Note that in both the 
cases beyond the freezeout surface the fraction of e^ pairs in the wind is negligibly small 
[39]. 

Note that our calculations probably do not apply to quark stars in the color flavor 
locked (CFL) superconducting phase. Previously it was suggested [40] that, despite the 
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absence of electrons in the bulk SQM in the CFL phase, the electrosphere may exist due 
to the surface quark charge [41]. However, the recent analysis [42] gives evidence in favor 
of absence of such a surface charge. But for the CFL phase may exist a significant photon 
emission from the SQM itself due to the photon-gluon mixing [43]. The results of [43] 
show that this radiation is comparable to the black body limit. Since we also obtain the 
radiation rate comparable to the black body radiation it may be difficult to distinguish a 
bare quark star in the CFL phase from that in normal (or 2SC) phase. 

5 Conclusion 

In summary, we have evaluated the photon emission from the electrosphere of a bare quark 
star (in normal or 2SC phase). The analysis is based on the LCPI reformulation [17] of 
the AMY formalism [16] to the photon emission from relativistic plasmas. The developed 
approach, contrary to the previous qualitative studies [13, 14, 15], allows, for the first time, 
to give a robust treatment of the Pauli-blo eking effects in the photon bremsstrahlung. We 
demonstrate that for the temperatures T ~ 0.1 ^ 1 MeV the dominating contribution to 
the photon emission is due to bending of electron trajectories in the mean electric field 
of the electrosphere. The energy fiux from the mean field photon emission is of order 
of the black body limit. Our results show that the contribution of the Bethe-Heitler 
bremsstrahlung due to electron-electron interaction is negligible as compared to the mean 
field photon emission. In contrast with [13] we demonstrate that the LPM suppression is 
negligible. 

The energy fiux related to the mean field bremsstrahlung turns out to be larger than 
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that from the tunnel e^e^ pair creation [4, 10] as well. In the light of these results 
the situation with distinguishing bare quark stars made of the SQM in normal (or 2SC) 
phase from neutron stars may be more optimistic than in the scenario with the tunnel 
e"'"e~ creation discussed in [25]. 
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Appendix A. Calculation of the function D{q_j_) 

In this appendix, we discuss the calculation of the function Z}(q_|_). To evaluate this 
function one needs to know the correlator D'^". In momentum representation one can 
obtain 

D^'^(g) = -2[l + nB(go)]ImDr(?), 

where ns = [exp{qQ/T) — 1]~^ is the Bose-Einstein factor, and D!^'^{q) retarded Green's 
function. As was already noted the function P{p) is gauge invariant, and one can use V^" 
in any gauge. Expressing the retarded propagator in the Coulomb gauge in terms terms 
of longitudinal and transverse photon self-energies one can obtain 



^. ^ = _i 7 a exp(go/T) [ Imn^(go,q) 

^"^^^ n J ^°exp(an/T) - 1 [q^ - ReH, fan, q)!^ + (h 



'exp(go/T) - 1 [ [q^ - ReH^lgo, q)]^ + (Imn^lgo, q))^ 

qi imnT(go,q) 



+ 



q2 [qi + RenT(go, q)]' + (ImnT(go, q))' 
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(A.l) 



qz=qo 



In numerical calculations we use for IIl,t the HDL expressions [30, 31] 



nL(go,q) 



m 



D 



^In *^ -1 

2g \qo-qJ 



nT(go,q) 



m 



D 



% , {T-%)%^ (Qo + q 



2g3 



In 



.qo-q, 



(A.2) 
(A.3) 



with the Debye mass m|, = ^(/i^ + -k^T^/3). 



Appendix B. Formulas for the light-cone wave func- 
tions 

For zero f the light-cone wave functions have definite orbital quantum number m. As 
was mentioned the light-cone wave functions appear from the longitudinal integrals of the 
Green's function. For f = it is the free Green's function given by 



^1.(6, P2l6> Ply 



M 

2mi 



exp < I 



M{p, - p,Y ie' 



2C 



2M 



(B.l) 



with C = ^2 ^ Ci- The ^-integration can be performed with the help of the relation 



iM 

rf^/C^(p2,0|pi,O = ^o(|p2-Pi|e), 

-oo TT 



(B.2) 



where Kq is the Bessel function. Then the light-cone wave functions can be written in 
terms of the Bessel functions Kq and Ki. After representing the vertex operator (6) in 
terms of the helicity state projectors as in [17] one can obtain for Ae' = Ag 



1 fW 
^(a;,p, Ae,Ae',A^) = --W-- [A^(2 - x) + 2Aea;] exp(-iA^v9)ei^i(pe) , 

ZTT V IX 



where (/? is the azimuthal angle. For Ae' = — Ag 



(B.3) 



—I 



^(x, p, Ae, -Ae, 2Ae) = -—V2ax^meKo{pe) . 

Zn 



(B.4) 
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Figure 1: The function P(p) (9) in units of the Debye mass versus prnD for different values 
of the ratio r = T/mo- (a) shows the total L + T contribution, (b) and (c) show L and T 
contributions, respectively. The curves correspond to: r = (solid line), r = 0.5 (dotted 
line), r = 1 (short dashes), r = 2 (long dashes). The thick solid line in panel (a) shows 
prediction of the static model obtained with the dipole cross section (11). 
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Figure 2: The contributions to the spectrum dP/dxdL from the mean field mechanism 
(solid line) and the fiuctuation mechanism (dashes) for /i = 10 MeV at T = 0.2 and 
T = 1 MeV. The thick curves are for nonzero photon mass, and the thin ones are for 
massless photon. The contribution of the fiuctuation mechanism is calculated using the 
Bethe-Heitler term with the distribution (24). 
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Figure 3: The emissitivity versus the photon energy u for fi = 5 and /i = 10 MeV at 
T = 0.2 and T = 1 MeV. The thick sohd hne shows the mean field bremsstrahlung. The 
contribution of the fiuctuation mechanism is shown without (dashes) and with (thin solid 
line) the mean field suppression. 
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Figure 4: The differential radiated energy flux from the electrosphere for the mean fleld 
bremsstrahlung (thick sohd hne) and for the Bethe-Heitler bremsstrahlung with (thin 
sohd hne) and without (dashes) the mean fleld suppression. The dotted curves show the 
black body spectrum. 
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Figure 5: The total radiated energy flux (scaled to the black body radiation) from the 
electrosphere for the mean fleld bremsstrahlung (thick solid line) and for the Bethe-Heitler 
bremsstrahlung with (thin solid line) and without (short dashes) the mean field suppres- 
sion. The contribution from the the tunnel e~^e~ creation [4, 10] evaluated using (25) is 
also shown (dotted line). The long dashes show the results for e^ + e^ — )■ e^ + e^ + 7 
process obtained in [13]. The dot-dashed line show the results for the same process of 
[14]. The dot-dot-dashed line shows the bremsstrahlung contribution with inclusion of 
the mean Coulomb field of [15]. qc 



